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Abstract. Studies of the molecular interstellar medium that fuels star formation and 
supermassive black hole growth in galaxies at cosmological distances have undergone 
tremendous progress over the past few years. Based on the detection of molecular gas 
in >120 galaxies at z=l to 6.4, we have obtained detailed insight on how the amount 
and physical properties of this material in a galaxy are connected to its current star 
formation rate over a range of galaxy populations. Studies of the gas dynamics and 
morphology at high spatial resolution allow us to distinguish between gas-rich merg- 
ers in different stages along the "merger sequence" and disk galaxies. Observations 
of the most massive gas-rich starburst galaxies out to z>5 provide insight into the role 
of cosmic environment for the early growth of present-day massive spheroidal galax- 
ies. Large-area submillimeter surveys have revealed a rare population of extremely 
far-infrared-luminous gas-rich high-redshift objects, which is dominated by strongly 
lensed, massive starburst galaxies. These discoveries have greatly improved our un- 
derstanding of the role of molecular gas in the evolution of massive galaxies through 
cosmic time. 



1. Introduction 

Great progress has been made in studies of galaxy evolution out to high redshift over the 
past years, but there are a number of fundamental questions that remain to be answered. 
One of the most important remaining issues is to understand whether star formation 
and the subsequent buildup of stellar mass in galaxies at early cosmic times occurs 
dominantly through major mergers (e.g., Springel et al. 2005) or through a combination 
of minor mergers and steady, so-called "cold-mode" accretion (e.g., Dekel et al. 2009). 
This issue is closely connected to the question how galaxies obtain the (dominantly 
molecular) gas that fuels star formation, and what their gas mass fractions are. Gas 
mass fractions, in relation to star formation rates, determine the evolutionary state of a 
galaxy (e.g., Daddi et al. 2010a; Tacconi et al. 2010). Next to these more general issues, 
it is important to better understand the physical properties, chemical composition, and 
dynamics of the star-forming gas in high-z galaxies, which set the initial conditions for 
star formation (e.g., Riechers et al. 2006a, 2008a; WeiB et al. 2007). 

All of these fundamental issues are directly tied to studies of molecular gas in high 
redshift galaxies. Despite the great progress that has been made in this field over the 
past 15 years, our understanding was ultimately limited by the technical capabilities 
of past observatories in the centimeter to submillimeter wavelength range, where the 
most common molecular gas diagnostics can be observed at high redshift. New ob- 
servatories, such as the Atacama Large sub/Millimeter Array (ALMA) that currently 
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nears completion, will be key to ultimately solve many of the remaining mysteries in 
this field. This article summarizes some of the most recent (pre- ALMA) progress in the 
field of molecular gas observations out to high redshift, and how the gas properties en- 
tangle with the driving mechanism of star formation in massive high redshift galaxies. 

2. CO Detections at High Redshift: a Brief Summary 
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Figure 1 . Detections of CO emission in z> 1 galaxies. Left: Total number of detec- 
tions (red) and detections per year (yellow) since the initial detection in 1991/1992. 
Right: Detections as of December 201 1 as a function of redshift, and color encoded 
by galaxy type (figure updated from Riechers 201 la). 

To date, molecular gas (most commonly CO) has been detected in >120 galaxies 
at z>l (see reviews by Solomon & Vanden Bout 2005; Omont 2007; Fig. 1), back to 
only 870 million years after the Big Bang (corresponding to z=6.42; e.g., Walter et al. 
2003, 2004; Bertoldi et al. 2003; Riechers et al. 2009). The bulk of these galaxies are 
massive, hosting molecular gas reservoirs of Mgas>10^'' M©, commonly with high gas 
fractions of (at least) tens of per cent. The sensitivity of past observatories has only 
allowed us to probe less massive and/or less gas-rich systems with the aid of strong 
gravitational lensing (e.g.. Baker et al. 2004; Coppin et al. 2007; Riechers et al. 2010a, 
2011a; Swinbank et al. 2010, 2011). 

Approximately 20% of the detected systems are massive, gas-rich optically/near- 
infrared selected star forming galaxies (SFGs; e.g., Daddi et al. 2010a; Tacconi et 
al. 2010), and 30% each are far-infrared-luminous, star-bursting quasars (QSOs; e.g., 
Wang et al. 2010; Riechers 201 lb) and submillimeter galaxies (SMGs; e.g., Greve et al. 
2005; Tacconi et al. 2008; Fig. 1). The rest of CO-detected high-redshift galaxies are 
limited samples of galaxies selected through a variety of techniques, such as Extremely 
Red Objects (EROs), Star-Forming Radio-selected Galaxies (SFRGs), 24 jum-selected 
galaxies, gravitationally lensed Lyman-break galaxies (LBGs), and radio galaxies (RGs; 
see Riechers 201 la for a recent summary). 

Besides CO, the high-density gas tracers HCN, HCO"^, HNC, CN, and H2O were 
detected towards a small subsample of these galaxies (e.g., Solomon et al. 2003; Riech- 
ers et al. 2006b, 2007, 2010b, 2011b; Guelin et al. 2007; Omont et al. 2011). 

3. Luminosity Relations and the Star Formation Law 

The CO luminosity L^q is commonly considered to be a measure for the total molecular 
gas mass Mgas in galaxies, and the far-infrared (FIR) luminosity LpiR is considered to 
be a measure for the star formation rate (SFR; e.g., Solomon & Vanden Bout 2005). 
Thus, the relation between L^q or Mgas and Lfir or SFR may be considered a spatially 
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integrated version of the Schmidt-Kennicutt "star formation law" between gas surface 
density and star formation rate surface density (e.g., Kennicutt 1998). 




Figure 2. Comparison of CO luminosity/gas mass with (far-) infrared luminosity 
as a tracer of the star formation rate in low- and high-z galaxies. Left: Mgas-LiR re- 
lation for low-z spirals, starbursts and ULIRGs, and high-z SFGs ("BzKs/normal"), 
SMGs, and quasars. The solid line represents a fit to spiral galaxies and SFGs, and 
the dotted line shows the same trend offset by 1 . Idex in Lir (which is consistent with 
the starburst galaxies, ULIRGs, SMGs and quasars). These two lines may represent 
two sequences for disk and starburst galaxies, with the offset being due to different 
dynamical timescales for star formation (e.g., Daddi et al. 2010b). Right: L^q-Lfir 
relation for quasars at low and high z (fitted by the solid line), and galaxies without 
dominant AGN at low z (spirals and (U)LIRGs) and at high z (SMGs; fitted by the 
dashed line). At high L^q, galaxies with and without dominant AGN statistically 
occupy the same region, while at low L^q, there is tentative evidence for an off- 
set of quasars toward higher Lfir- This may indicate that AGN heating contributes 
significantly to Lfir at low L^q, but not at high L'^^ (Riechers 201 lb). 

Even in galaxies with luminous active galactic nuclei (AGN), LnR is commonly 
used as a proxy for the star formation rate in the host galaxy. In principle, both the 
AGN and star formation can heat the dust that gives rise to the continuum flux observed 
in the FIR, but the characteristic dust temperatures of AGN heating are typically by a 
factor of a few higher than those of dust heated by young stars. Thus, the warm dust 
in AGN-starburst systems observed in the rest-frame FIR is commonly thought to be 
dominated by heating within the host galaxy, in particular in the most intense, dust- 
enshrouded starbursts. If, however, Lfir were to be dominated by the AGN, one would 
expect an elevated Lfir for such galaxies in the L^^-Lfir relation. In Figure 2 (right 
panel), a comparison of the L^q-Lfir relation for nearby and high-z quasars to nearby 
galaxies, (ultra-)luminous infrared galaxies ((U)LIRGs) and SMGs without dominant 
AGN is shown (Riechers 2011b). For galaxies with low L^q, there is an indication for 
an excess in Lfir for quasars relative to other systems; however, there is no evidence 
for such a trend at high L^^. This may suggest that, in systems with relatively low gas 
and dust content, AGN contribute significantly to Lfir, but not in systems with high gas 
and dust content. Thus, for massive high-redshift galaxies, Lfir appears to be a good 
proxy for the SFR, even in quasars (Riechers 2011b). 

Beyond the issue of dust heating, the question occurs if star formation progresses 
the same way in all types of galaxies. In disk-like spiral galaxies like the Milky 
Way, star formation takes place in molecular clouds with compact, dense cores, con- 
fined by self gravity (e.g., Solomon et al. 1987). In mergers, such as the Antennae 
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(NGC 4038/39; Fig. 3), the gas and star formation appear to peak on the galaxy nu- 
clei, but also on relatively large scales in the dense overlap region between the merging 
galaxies, leading to the formation of so-called super star clusters (e.g., Wilson et al. 
2003). Such constellations are also commonly found in the highest-resolution CO stud- 
ies of high-redshift FIR-luminous massive galaxies (Fig. 3; e.g., Riechers et al. 2008a, 
2008b). In the nuclei of ULIRGs, the most extreme nearby starbursts, star formation 
appears to occur in a dense, intercloud medium, bound by the potential of the galaxy, 
rather than in virialized clouds (e.g., Downes & Solomon 1998). The differences be- 
tween starburst and disk galaxies are reflected in the star formation law. When compar- 
ing the Mgas-Z^iR relation for three largest CO-detected samples at high z, i.e., quasars, 
SMGs, and SFGs, to low redshift galaxies, two interesting trends occur. First, SFGs 
extend the relation found for nearby spiral galaxies to higher Mgas. Second, quasars 
and SMGs extend the trend found for the most intense nearby starbursts and ULIRGs 
to higher Mgas. Both trends agree with the same slope, but are offset by l.ldex in Ljr. 
Daddi et al. (2010b) interpret this as evidence for two sequences in this relation for disk 
galaxies and starbursts, which emerge from the different dynamical timescales of star 
formation in these systems (Fig. 2, left; see also Genzel et al. 2010). Recent theoret- 
ical studies have attempted to understand the different trends based on the underlying 
conversion factor, aco, from L^q to Mgas, and find that the two sequences could only 
be unified to a single relation when assuming a broad continuum of conversion factors 
(e.g., Narayanan 2012, this volume). 
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Figure 3. Molecular gas distribution in major mergers. Left: The z=4.41 AGN- 
starburst system BRI 1335-0417 (0.16"x0.14" resolution, corresponding to l.Okpc; 
Riechers et al. 2008a). Right: The nearby major merger NGC 4038/39 (the Antennae; 
convolved to l.Okpc resolution and rotated for illustration; Wilson et al. 2003). Both 
systems are dominated by a clumpy, asymmetric gas distribution, and some diffuse 
emission is seen at low flux density level. It is plausible to assume that the clumps 
can be identified with two galaxy nuclei and a dense overlap region in both cases. 

4. Gas Dynamics: A "Merger Sequence" of High Redshift Galaxies 

The advent of the Expanded Very Large Array (EVLAfl has recently enabled studies of 
the full gas content, distribution and dynamics of different populations of high-redshift 
galaxies through observations of €0(7=1— >0) emission (e.g., Riechers et al. 2010a, 
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2011c, 201 Id, 201 le; Ivison et al. 2011). In contrast to previous studies at shorter 
wavelengths in higher rotational lines of CO, observations of CO 7=1— >0 trace the full 
amount and extent of the molecular gas, and can be more directly compared to obser- 
vations in the nearby universe, which are most commonly undertaken in CO 7=1^0. 

One particular aspect of these studies is that they provide deeper insight into the 
mechanisms that are driving the conversion of gas into stellar mass in the most intensely 
star-forming galaxies at early cosmic times, such as SMGs. It has been argued in the 
past that most SMGs are major mergers of two gas-rich galaxies (e.g., Tacconi et al. 
2008; Engel et al. 2010). Detailed studies of CO(7=1^0) emission with the EVLA 
confirm this picture, and show a range in merger properties and stages (e.g., Riechers et 
al. 2011c, 201 Id; Fig. 4). Indeed, it becomes possible to place SMGs along a "merger 
sequence" of high redshift galaxies, ranging from systems with two disk-like gas-rich 
galaxies that are separated by tens of kiloparsec in projection and several hundreds 
kms~^ in velocity over actively merging systems with a single, morphologically and 
dynamically complex gas distribution to systems that have almost reached coalescence. 
These early investigations demonstrate that studies of larger samples at higher spatial 
resolution with the fully upgraded EVLA in the future will allow us to distinguish be- 
tween different galaxy populations based on the physical properties of their interstellar 
media. This will yield a more complete understanding of the processes that trigger star 
formation and black hole activity in the early universe, and their relative importance for 
the buildup of stellar mass in galaxies as seen at present day. 
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Figure 4. CO(7=1^0) observations of submillimeter galaxies along the "merger 
sequence" with the Expanded Very Large Array (Riechers et al. 2011c, 2011d). 
SMGs show complex gas morphologies that can extend over >10kpc scales, and 
commonly consist of multiple components. The gas distribution and kinematics of 
the majority of SMGs are as expected for major, gas-rich mergers. The observed 
diversity in these properties is consistent with different merger stages, as shown here 
for three examples. Left: The two gas-rich galaxies in this SMG are separated by 
tens of kpc and ~700kms"', representing an early merger stage. Middle: This SMG 
still shows two separated components, but at similar velocity, representing a more 
advanced merger stage. Right: This SMG shows a single, complex extended gas 
structure with multiple velocity components, representing a fairly late merger stage. 



5. Cosmic Environments of Massive Starbursts at Very High Redshift 



Recent studies have led to the discovery of the long sought-after high-redshift tail of 
SMGs at z>4 (e.g., Capak et al. 2008; Daddi et al. 2009). SMGs at the highest redshifts 
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are of particular interest, as they trace some of the most massive and active systems at 
early cosmic times. In cosmological simulations, such systems are expected to grow 
in the rare, most overdense regions in the early universe (e.g., Springel et al. 2005). 
Thus, very high-z SMGs may trace the most distant proto-cluster regions in the uni- 
verse, which have the potential to grow into the most massive cosmic structures seen 
at present day. Indeed, following the discovery of a SMG at an unprecedented redshift 
of z=5.298, observational evidence has recently been found for this mode of structure 
formation in the early universe (Capak et al. 201 1; Riechers et al. 2010c; Fig. 5). Deep 
imaging and spectroscopy of its cosmic environment as part of the Cosmic Evolution 
Survey (COSMOS) reveals a high overdensity of Lyman-break galaxies at similar red- 
shifts within a co-moving volume of 2 Mpc radius around the SMG. The structure of 
the central proto-cluster encompasses a halo mass of >4x1O''M0, which would be 
consistent with what is expected for the early formation of a present-day galaxy cluster. 
Clearly, future discovery and studies of SMGs at comparable and higher redshifts are 
desirable to determine the fraction of the first massive starbursts that are associated with 
such galaxy overdensities, as necessary to study the effects of environment on galaxy 
formation within a more general cosmological context. 
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Figure 5. Gas properties and environment of the most distant SMG AzTEC-3 
at z=5.298 (Riechers et al. 2010c; Capak et al. 2011). Left: 2'x2' region around 
AzTEC-3. The white circles show photometrically identified candidate Lyman-break 
galaxies (LBGs) at z=5.3, and the green circle indicates a region of 2 Mpc co-moving 
radius. The labels indicate spectroscopic redshifts of confirmed LBGs or of low- 
z galaxies and stars rejected based on their spectral energy distributions. Middle: 
Zoom-in of the central 0.865 Mpc (co-moving) region, with spectroscopic redshifts 
labeled. Right: Zoom-in of the SMG, with CO(y=5— >4) contours overlaid (top), 
and optical spectra of the SMG and companion galaxies (bottom). The massive, 
very high z starburst galaxy AzTEC-3 traces a "proto-cluster" region of high galaxy 
overdensity in the early universe, showing that the highest-redshift SMGs may be an 
ideal tracer of the most massive cosmic sttuctures at early cosmic times. 



6. Herschel: The Lensing Revolution 



One of the most intriguing recent discoveries made by the Herschel Space Observatory 
is the identification of a rare population of extremely luminous high-redshift starburst 
galaxies at submillimeter wavelengths (550()/jm>100mJy; e.g., Negrello et al. 2010). 
This population far exceeds the expected number counts of SMGs at the bright end, and 
is dominated by gravitationally lensed SMGs. This discovery is interesting, because it 
allows for a very efficient selection of gravitational lenses based on a flux density mea- 
surement alone (after rejection of low-redshift contaminants). Also, the natural mag- 
nification in size and flux density provided by gravitational lensing have enabled new 
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techniques to better understand the SMG population itself. First, the high expected CO 
emission line fluxes, paired with recent spectral bandwidth upgrades of radio/millimeter 
observatories, have yielded a significant number of "blind" CO detections in these and 
similar sources without any prior constraints on their redshifts (e.g., Riechers 2011b; 
Frayer et al. 2011; Cox et al. 2011; Harris et al. 2012; Lupu et al. 2012). Second, high- 
resolution CO imaging studies, in combination with detailed lens modeling, allow us to 
probe down to structure sizes that would remain inaccessible with current-generation 
instruments without the aid of gravitational lensing (e.g.. Fig. 6; Riechers et al. 2011a; 
Swinbank et al. 2011; see also Riechers et al. 2008a). Only ALMA will match the 
sensitivity and resolution of these studies in unlensed galaxies, and will even allow us 
to probe down to individual molecular cloud scales in lensed systems. 




Figure 6. Multi-wavelength observations and modeling of HLSW-01, a rare, 
exceptionally bright, lensed SMG at z-2.951 discovered in the HerMES Her- 
ic/ieZ/SPIRE survey (Riechers et al. 2011a; Gavazzi et al. 2011). Left: Overlay of 
CO(7=5— »4) contours on observed-frame 890 ^m. continuum emission. A to D label 
the lensed images of the galaxy. Middle: Same overlaid on a 2.2 /im image. Right: 
Lens model of the 2.2 /vm {left) and CO(7=5— >4) emission {right) in the image {top) 
and source {bottom) planes. This SMG is lensed by a galaxy group at z^O.6, yielding 
an extreme apparent 250 ;um flux density of 425+lOmJy due to a lensing magnifica- 
tion factor of ~10.9, and a wide lens image separation of ~9" . This source represents 
the bright end of a new population of lensed SMGs discovered with Herschel. 

7. Summary and Outlook 

Studies of the molecular interstellar medium in galaxies at cosmological distances have 
shown great progress in the past few years. These advances were possible due to a 
combination of improved selection techniques to identify gas-rich galaxies in the early 
universe and major improvements in instrumentation, which have enabled studies of 
less extreme galaxy populations than previously possible, and have liberated high red- 
shift molecular line studies from the prerequisite of a precise spectroscopic redshift ob- 
tained through other measures (which were a major selection bias in the past). The most 
recent studies thus finally probe beyond the "tip of the iceberg" of the most far-infrared- 
luminous galaxies with optical spectroscopic redshifts that were the focus of molecular 
line studies at high z in the past. These observations also offer a "sneak peek" into the 
detailed investigations of the physical properties of galaxies at early cosmic times that 
will become possible with ALMA, once it commences full science operations. 
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